Radiation reaction effects on ion acceleration in laser foil interaction 
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Radiation reaction effects on ion acceleration in laser foil interaction are investigated via ana- 
lytical modeling and multi-dimensional particle-in-cell simulations. We find the radiation effects 
are important in the area where some electrons move backwards due to static charge separation 
field at the laser intensity of 10^^ W/cm^. Radiation reaction tends to impede these backwards 
motion. In the optical transparency region ion acceleration is enhanced when the radiation effects 
are considered. 
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INTRODUCTION 

Along with the fast development of laser technology, 
the power intensity which is inaccessible one decade be- 
fore becomes reality. The focused laser intensity of 
10^^ W/cm? is available today. The new generation 
of laser system such as ELI is under construction. Fo- 
cused intensity of 10^^ W/cm? is just forward jl|. By 
use of such intense laser power, plenty of physical prob- 
lems can be studied and amount of applications are wait- 
ing for discovering. Among them particle acceleration 
presently and will still attract extensive attentions. Elec- 
tron acceleration by laser plasma accelerator is aimed to 
TeV [2I, ion acceleration also towards to stable, high and 
mono-energetic beam [3]. Besides ultrahigh power laser 
system, plasma mirror technology also makes the high 
contrast clean pulse possible [4], which makes ultrathin 
target (~nm thickness) now applicable for ion accelera- 
tion Unlike the target normal sheath field acceler- 
ation (TNSA) [6], radiation pressure dominated ion ac- 
celeration (RPA) can have a much longer acceleration 
distance [7]. Simulation results show ions can get GeV 
energy when a 20fs long laser pulse with 10^^ W/cm? 
intensity is used [8]. Besides the two mentioned ion ac- 
celeration mechanisms, recently ion acceleration in the 
optical transparency region named the "laser breakout 
afterburner" (BOA) ji| has also been studied. It demon- 
strates that ions can get continuously acceleration even 
the laser pulse transmits through the ultrathin foil, which 
also gives a possible way to get high energy ion beams. 

As is well known, high energy particles when suffer- 
ing acceleration self-radiation is concomitant. Radiation 
reaction effects should be considered when the radiation 
damping force is comparable to the external one. For the 
above mentioned processes for particle acceleration, elec- 
trons are always endured intensive acceleration. When 
the laser intensity increases further, the radiation reac- 
tions cannot be neglected. Naumova et al have stud- 
ied the radiation reaction effect on the laser hole boring 
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Figure 1: (color online) Sketch map of laser electron interac- 
tion. Here a linearly polarized laser pulse is used. 



process and point out it plays a positive role as it al- 
lows one to maintain the electron thermal energy on a 
relatively low level and prevents the electron backward 
motion through the pulse (iq| . 

In the present paper, we study their effects on ion ac- 
celeration. As a primary result, we find the radiation 
mainly comes from the electrons moving backwards in 
the laser pulse and the radiation damping impedes this 
kind of backward motion, which can reduce the particles' 
volume in the phase space and improve the ion accelera- 
tion energy and quality. 



RADIATION AND ITS EFFECT 

Before showing the effect of radiation damping, we first 
check the threshold of the laser intensity for the impor- 
tant radiation reaction effects from existed fundamental 
formulae. For an electron with velocity v^c, one can get 
the radiation power P{t) at the local radiation time t as: 
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when dacce-^^' Here p is the electron momentum, 7 is the 
relativistic factor, e and m are the charge and mass of the 
electron, respectively. So the radiation due to the trans- 
verse acceleration is 7^ stronger than the one due to the 
longitudinal acceleration. For the synchrotron radiation, 
the radiation is mainly in the direction of electron motion 
and concentrated within an angle of 6 with A6 ~ -. 

Based on the above knowledge, in the PIC code we 
only consider the radiation due to the transverse accel- 
eration and only the electrons whose 7 > 5 are assumed 
to contribute the radiation. 

By use of the normalized variables as in the PIC code 
{p^p/mc^t^t/To — cjoV^tt), we can get the momentum 
variation rate as: 
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where h is the Planck constant, coq is the laser frequency 
and Epho is the photon energy of the laser pulse. To be 
simply, in the PIC code we also assume the radiation is 
in form of photons. The photon frequency is ujc and the 
corresponding photon number is Uph = Sne'^dp/dhc. The 
radiation reaction force on the electron then is: 
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The reaction force points to the opposite direction of the 
electron motion. On the other hand, the electron also 
feels the external force as (Here we only show the trans- 
verse force due to laser field.): 
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where E and B are the intensities of electric and magnetic 
fields normalized by muoQc/e. We can get the threshold 
of the laser intensity for the obvious radiation damping 
effect by using: 
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It is: 
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Usually in a laser pulse, |^| = \B\ = a. For an electron 
moving along with a linearly polarized laser pulse, we 
have: 

Py = a 
Pa, = a^/2>0 
7 = l + aV2 



Here a is the normalized laser field. To get an obvious 
radiation effect, it should satisfy 7a>8.47 x 10'', corre- 
spondingly the laser intensity should satisfy: a>550. It 
is about 4.2x10^^ W/cm'^ for a laser with a wavelength 
of 1 /im, which is higher than the current running laser 
system. However, for an electron with longitudinal veloc- 
ity of 0, Eq. jSl) changes to j'^a > 8.47 x 10^. When 
the electron moves in the opposite direction of the laser 
pulse with the longitudinal velocity /3x— — 1, then one 
can get 7^a > 4.23 x 10^. That is the threshold of a can 
be 7 times smaller when the electrons with the same en- 
ergy (7) move backward (see Fig. [T]). The reason is that 
the radiation force is proportion to the square of the ac- 
celeration force. For the electrons moving forwards the 
transverse force they feel is about zero {E -\- P^^B ^ 0), 
so the radiation force is even smaller. However, for the 
electrons moving backwards, the felt transverse accelera- 
tion force is about AttQcE, so the radiation force is much 
larger. 



SIMULATION RESULTS AND DISCUSSION 

In the following, we use PIC simulations to check the 
radiation reaction effects and show their importance on 
ion acceleration. To simulate the damping effects we 
suppose that, at any given moment of time, the elec- 
tron radiation spectrum is synchrotron like (lll|. The 
critical frequency coc is given by the relation coc = 
(3/47r)7^|APx|/((it); AP± is the variation of transverse 
electron momentum force during the time step of dt. In 
our PIC code, we follow trajectories of each electron and 
calculate the emission during the interaction. We cal- 
culate the damping effects by considering the electron's 
recoil due to the emitted radiation. The recoil force is 
included in the equations of electron motion. We should 
mention that our method is different with the one used 
by Martins et al in OSIRIS code, in which they can 
get the radiation field and frequency in a faraway vir- 
tual detector pJi]. In our code we do not pay attention 
to the received radiation field on the virtual detector. 
Our main attention is focused on the reaction effects on 
plasma itself. The radiation we recorded every simulation 
step is the photons radiated in the local time of radia- 
tion, not the one at the observation time. Our method 
is also simpler than the one used by Sokolov et al, in 
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Figure 2: (color online) (a) Spatial-temporal distribution 
of the electron density, (b) Electron energy spectra in the 
simulations with and without radiation damping effects at 
t = 20 To. (c) The energy angular distribution of the radi- 
ated photons during t = 20 To and t = 21 To. The colorbar 
shows the relative photon number in the logarithmic value, 
(d) Proton spectra in the simulations with and without ra- 
diation damping effects at t = 20 To. Here a = 100 and 
n = 100. 
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Figure 3: (color online) (a) Spatial-temporal distribution 
of the electron density for the simulation without radiation 
damping effects, (b) The case with radiation damping effect, 
(c) The energy angular distribution of the radiated photons 
during t = 20 To and t = 21 To. The colorbar shows the 
relative photon number in the logarithmic value, (d) The 
radiated photon energy and positions during t = 20 To and 
t = 21 To. The colorbar shows the relative photon number in 
the logarithmic value. Here a = 180 and n = 100. 



which the modified non-perturbative Lorentz-Abraham- 
Dirac equation is resolved for particle motion instead of 
the normal Lorentz equation Although with our 

method the code cannot gives a correct radiation field 
on a virtual detector, it is relatively simple and can give 
appropriate description for the plasmas under radiation 
damping. Our findings are similar with Naumova et al. 
as they study the hole boring process for fast ignition, 
however, our main interest is on the ion acceleration in 
the laser foil interaction. 

We take ID-PIC simulations with the KLAP-ID 
code [l3|. In the simulation, the laser pulse has 
trapezoidal temporal intensity profile ("linear increase- 
plateau-linear decrease") with the pulse length lA/c- 
16A/c-lA/c. Here A = 1 /im is the laser wavelength. 
The normalized laser electric field a = eE/mcooc chang- 
ing from 100 to 500 are used. The target plasma are 
composed of electrons and protons and they are initially 
uniformly distributed from x = 4 A to x = 4.3 A. The 
plasma density is fixed to be 100 ric with the critical den- 
sity Tie = l.lxlO^^/cm^. 

For the presently fixed plasma density we find when 
the laser intensity is lower than a~125 the radiation ef- 
fects on the particle energy spectrum is not observable. 
Fig. [2Ka) shows the spatial temporal distribution of the 
electron density in the simulation of a = 100. In this 
condition, the ions are accelerated in the radiation pres- 
sure dominated region. Electrons and ions are moving 
together and the acceleration is phase stable. No obvious 
electron backward motion happens. The electron energy 



spectrum does not vary so much between the two simu- 
lations with and without radiation damping as shown in 
Fig-EJc). There is a sharp peak in the proton spectrum 
with the position at 0.63 GeV after 16 To acceleration 
time which is close to the theoretical value of 0.66 GeV. 
Fig. [2](c) shows the energy angular distribution of the ra- 
diated photons. Most of the radiation is in the forward di- 
rection and the photon energies are about 10 KeV. From 
the simulation we calculate the average radiation power 
is about 1.64 x 10^'' W/cm'^ which is negligible compared 
to the laser power of 2.76 x 10^^ W/cm'^. 

When we increase the laser intensity further to a = 
180, the radiation damping effects obviously appears. 
Fig. [3Ka,b) show the spatial temporal distribution of the 
electron density in the two compared simulations (with 
and without radiation effects) . As we see when there is no 
radiation effects amount of electrons move backwards and 
the electron density distribution disperses in space. On 
the contrary, after considering the radiation effects this 
kind of backward acceleration has been suppressed. More 
electrons are concentrated and move along with the ion 
bunch as shown in Fig. [3Kb). We show the energy angular 
distribution of the radiated photons in Fig.[3Kc). It com- 
poses of two parts. Except the lower part with low pho- 
ton energy which is similar with the one in Fig. ^c) (We 
call this "low energy radiation".), there is a high energy 
part whose radiation direction is uniformly distributed in 
the forward and backward directions (We call this "high 
energy radiation"). To find the source electrons for these 
radiated photons, we show the radiation position distri- 
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Figure 4: (color online) (a) Electron distribution in the x—px 
phase space, the momentum is normalized by rrieC] (b) Proton 
distribution in the x — px phase space; (c) Electron energy 
spectrum, the momentum is normalized by m^c; (d) Proton 
energy spectrum. The time here is t = 20 Tq. 



bution in Fig.[3Kd). Comparing with the positions of the 
electrons at t = 20 To in Fig. [S^b), we can see the low 
energy radiation {Ephoton < 1 MeV) mainly comes from 
the electrons moving with the ion bunch and the radi- 
ation mainly concentrates within the angle of 20^, the 
high energy photons {Ephoton^'^'^'^ MeV) are radiated 
by the electrons which move behind the ion bunch where 
the plasma is transparent to the laser pulse and electrons 
move both forwards and backwards, the radiation is also 
almost uniformly in the angular distribution. The to- 
tal radiation power density is about 8.84 x 10^^ W/cm'^ 
which is almost 10% of the laser power density. 

Particle distributions in the phase space are shown in 
Fig. [4Ka,b). Electrons have much smaller volume in the 
phase space when the radiation damping is included in 
the simulation. In this case no backward ion acceleration 
has been found. Much more ions are confined in the front 
accelerating bunch compared with the simulation with- 
out radiation damping. The electron energy spectrum 
also shows that electrons have a much lower peak energy 
(10 MeV) in the radiation included simulation, however, 
it is about 150 MeV in the simulation without radiation 
damping. Proton energy spectrum shown in Fig. [H con- 
firms much more protons are concentrated and acceler- 
ated in the bunch whose peak energy is about 1.5GeV at 
t = 20 To. 

From the energy point of view, at t = 20 To, to- 
tally 5.01 X 10^ J /cm? laser energy density has been 
transported into the simulation box. When no radia- 
tion damping is included, among them 2.44% (1.22 x 
10^ J I err?) transforms to electrons and 6.93% (3.47 x 
10^ J I cm?) transforms to protons. When radiation 
damping is included, 1.23% (6.15 x 10^ J /err?) trans- 



forms to electrons and 11.16% (5.59 x 10^ J / err?) trans- 
forms to protons. As we see radiation damping reduces 
the electron energy, however, it improves the proton ac- 
celeration. Neglect of radiation damping gives a lower 
estimation of energy conversion efficiency and worse spec- 
trum. 

In our simulations, we find that the high energy ra- 
diation appears when a > 125. Its quotient among the 
total radiation increases with the laser intensity. When 
a > 300, there is only high energy radiation. These dif- 
ferent radiation profiles also refiect the different accelera- 
tion scenarios. From the balance between the forces due 
to electrostatic field and radiation pressure 
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where / is the thickness of the target normalized by laser 
wave length (Ao) and /3e = a/{a -\- y/rriino) is the rela- 
tivistic hole boring velocity [15], we can get the critical 
laser intensity for the laser pulse transmitting through 
the target. For our present simulation parameters: no = 
100, / = 0.3, rrii = 1836, we get acr ^ 117. The intensity 
is close to the critical value (a ^ 125) for the high energy 
radiation. It means when the laser intensity is higher 
than acr J the target tends to be transparent to the laser 
pulse. Stable (at least in the ID case) structure of radia- 
tion dominated ion acceleration begins mixing with other 
acceleration mechanism (such as BOA). However, some 
ions can still be accelerated by the radiation pressure. 
The amount of these ions depends on the laser intensity. 
Other ions are accelerated or decelerated by the dispersed 
electrons. When the laser intensity increases further, ra- 
diation pressure acceleration disappears completely. Ions 
can only be accelerated in the heated electron pool which 
moves with the laser pulse. In this scenario, the radia- 
tion damping is important and necessary, which actually 
makes the electron pool cooling down and improves the 
ion acceleration. Further investigation still needs for this 
kind of acceleration scenario. 



SUMMARY AND DISCUSSION 

In summary, by use of PIC simulations we studied the 
radiation reaction effects on the ion acceleration in the 
radiation pressure dominated region and the transpar- 
ent plasma region. We find radiation damping effects are 
important for the electrons moving backwards and im- 
mersed in the laser pulse. Self-radiation impedes back- 
ward motion, cools down the electrons and makes more 
ions be concentrated and accelerated. We notice recently 
many studies on ion acceleration by use of ultrathin foil 
target, in which target plasma is almost transparent to 
the laser pulse. Electrons are dancing with the trans- 
mitted relatively long laser pulse. The local charge sep- 
aration field takes in charge of ion acceleration. In this 
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condition, a correct electron distribution in phase space 
is important to get the correct final maximum acceler- 
ated ion energy and acceleration scaling. When the laser 
intensity is larger than a = 100, radiation reaction could 
change the electron distribution in phase space. Further- 
more, with future laser system such as ELI (a > 1000), 
laser intensity is high enough to awake the contribution of 
radiation damping. Even for the electrons move forward, 
the radiation reaction should be considered. It deserves 
and is necessary to include the radiation reaction effects 
in the future PIC simulations when ultra intense laser 
pulse is used. 
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